Hot Quark Matter with an Axial Chemical Potential 
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We analyze the phase diagram of hot quark matter in presence of an axial chemical potential, 
pis. The latter is introduced to mimic the chirality transitions induced, in hot Quantum Chromody- 
namics, by the strong sphaleron configurations. In particular, we study the curvature of the critical 
line at small /x 5 , the effects of a finite quark mass and of a vector interaction. Moreover, we build 
the mixed phase at the first order phase transition line, and draw the phase diagram in the chiral 
density and temperature plane. We finally compute the full topological susceptibility in presence of 
a background of topological charge. 
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I. INTRODUCTION 

Understanding the properties of strongly interact- 
ing matter, in extreme conditions of high temperature 
and/or large baryon density, is very important to get a 
deeper knowledge of our universe, at the macroscopic as 
well as at the microscopic levels. For example, very hot 
matter, with estimated temperature of the order of 10 12 
Kelvin, is produced in heavy ion collisions (HICs) experi- 
ments immediately after the collision, see for example [l| 
for an indirect measurement of such a temperature by the 
PHENIX collaboration at Brookhaven National Labora- 
tory (BNL). As a consequence, it is crucial to make theo- 
retical investigations on the phase structure of Quantum 
Chromodynamics (QCD), the theory of strong interac- 
tions, in conditions of high temperature, in order to give 
proper interpretation of experimental data, at the same 
time suggesting new phenomena to look for. 

The theoretical knowledge on the thermodynamics of 
QCD at zero baryon density is fixed by Lattice simu- 
lations, which allow first principle numerical computa- 
tions of the relevant thermodynamical quantities. In fact, 
Lattice simulations performed by independent groups 
show that in the range of temperature (140, 180) MeV, 
a crossover from the hadron phase to the deconfincmcnt 
phase takes place, see 0, Q and references therein. On 
the other hand, at finite baryon chemical potential it is 
not possible to use first principle Lattice simulations to 
investigate the structure of QCD (with three colors), be- 
cause of the sign problem, see [1] for a review. Lattice 
simulations in the strong coupling limit seem promising 
to avoid the sign problem, see 0, 0, even if the prob- 
lem of the continuum limit within this approach has to 
be considered with care: in particular, the continuum 
limit is not reached in but the 1/g 2 evolution in a 
range of j3 = 6N c /g 2 < 4 = j3 s is consistent with Lattice 
Montecarlo simulations, for which the continuum limit is 
realized for j3 « 6, hence not too far from f3 s . 

The lack of feasible first principle calculations of the 
QCD thermodynamic properties at finite baryon chemi- 
cal potential invokes the use of other theoretical strate- 



gies to investigate the structure of QCD in this regime. 
The main strategy is the use of some model. Among the 
several models, the Nambu-Jona-Lasinio (NJL) model Q 
is very popular, see Q for reviews. In the NJL model, the 
QCD gluon-mediated interactions are replaced by effec- 
tive interactions among quarks, which are built in order 
to respect the global symmetries of QCD. Under some 
approximations, it is possible to derive the NJL model 
effective interaction kernel from first principles QCD, 
see 0, [T3| • The common feature of the effective models 
is that they share (most of) the symmetries of the QCD 
Lagrangian. At the same time, they are able to describe 
spontaneous chiral symmetry breaking, as well as other 
kinds of breaking patterns expected at large chemical po- 
tential, in a self-consistent way. Moreover, it has been 
shown how it is possible to extend the chiral models in or- 
der to compute quantities which are sensible to confinemt 
properties of a given phase [ill > see 0413, H3, [M| for re- 
cent studies. One of the advantage of these models is that 
they allow for analytic or semi-analytic computations at 
zero as well as at finite chemical potential, at least at 
the one-loop level; thus they do not rely (in general) on 
simulations, and are not affected by the complex value of 
the quark determinant, allowing speculations about the 
structure of QCD in regimes where first principles cal- 
culations are not feasible. After standard bosonization 
procedure, the chiral and deconfinement transitions are 
described in terms of collective fields which take some 
expectation value, and whose quantum fluctuations de- 
scribe physical particles (a and ir mesons in the simplest 
version of the model) . 

It has been suggested that very strong magnetic fields 
are produced during the very first moments of a noncen- 
tral heavy ion collision 2^- 27|; this has motivated sev- 
eral studies about the effect of a strong magnetic back- 
ground on the QCD phase structure, see [28H39j and ref- 
erences therein. Moreover, since the temperature of the 
fireball produced by the collision is very high, a copi- 
ous production of topological giuon configurations (i.e., 
the QCD sphalerons) with finite winding number is ex- 
pected [40], which induce locally chirality imbalance in 
the hot plasma as a natural consequence of the QCD 
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Ward identity. The combined effect of the latter and 
of the magnetic field induces an electric current along 
the direction of the magnetic field. This effect, called 
the Chiral Magnetic Effect (CME) (HEJ, leads to the 
event-by-event separation of electric charges with respect 
to the reaction plane, which is a parity {V) as well as a 
CP-odd effect. Experimental data obtained by the STAR 
collaboration at BNL seem to point towards the direction 
of charge separation in collisions [42j], even if the inter- 
pretation of such data in terms of the CME is still under 
debate [USES]- 

Because of the expected production of chirality imbal- 
ance in the quark-gluon-plasma phase of QCD, and be- 
cause of its potential relevance for the physics of heavy 
ion collisions, it is of interest to study how chirality modi- 
fies the structure of QCD itself. This is the main scope of 
this article. In particular, we continue the study of (45l — 
[48j in which chirality was induced by an axial chemi- 
cal potential, ^5, conjugated to chirality. Besides the 
applications to the CME, which are very interesting on 
their own because of the potential relevance for the phe- 
nomenology of heavy ion collisions, the theory at finite 
/Xg is interesting because it does not suffer from the sign 
problem; as a consequence, grandcanonical ensembles at 
finite fj>5 can be simulated on the Lattice [49|, l50( ■ These 
studies might be helpful to understand the structure of 
the QCD phase diagram at finite baryon chemical po- 
tential. For example, in [13] it has been suggested that 
the critical endpoint of the QCD phase diagram might be 
detected in Lattice simulations at finite ^5. The details 
about the theoretical framework will be given in Section 
II of the article. Here, we wish to stress the novelties 
embraced by our study. 

Firstly, we investigate the effect of the vector inter- 
action, as well as the finite current quark mass, on the 
location of the critical endpoint. Such aspects should be 
taken into account if a comparison with the Lattice data 
is desirable, but they have not been studied in [13]. More- 
over, it is of interest to analyze analytically the effect of 
the axial chemical potential on the chiral condensate, and 
on the curvature of the critical line for restoration of chi- 
ral symmetry. The latter aspect is quite interesting, since 
it shows how a competition between the vacuum term 
and the thermal excitations compete and eventually lead 
to a lowering of the critical temperature. This reduc- 
tion of the critical temperature was found numerically in 
the previous model studies [45l - [47j but it was not inves- 
tigated in detail. Thirdly, we are interested to the phase 
diagram in the 715 — T plane. This might be interesting 
because n§ is connected to the topological charge density, 
because of the integrated Ward identity. With an abuse 
of nomenclature, we call the latter as canonical ensemble 
formulation, in analogy with the case of QCD at finite 
baryon density. Finally, we compute the full topological 
susceptibility at zero as well as at finite as a function 
of temperature. The latter part of our investigation has 
some overlap with the model study of Ref. [19J, where 
topological susceptibility with a background of topolog- 



ical charge is computed. In the case of Ref. [T^, the 
background topological charge is introduced by adding a 
0— term to the QCD action. 

The plan of the paper is as follows. In Section II we 
introduce the axial chemical potential. In Section III 
we summarize the model we use in our calculations. In 
Section IV we perform a perturbative analysis at small 
/is at zero as well as at finite temperature, computing 
the dependence of the chiral condensate and of the criti- 
cal temperature on the axial chemical potential. In Sec- 
tion V we relate the phase diagrams in the canonical 
and grand-canonical ensembles, and build explicitly the 
mixed phase. In Section VI we study the effect of the 
bare quark mass, and of a vector interaction, on the lo- 
cation of the critical endpoint of the phase diagram in 
the 115 — T plane. Finally, we present our conclusions in 
Section VII. 

II. THE AXIAL CHEMICAL POTENTIAL 

In this Section we define the axial chemical poten- 
tial, /X5. It has been already discussed in several refer- 
ences, see [4l|, I45U47I l5lj and references therein. There- 
fore we limit ourselves to the basic definitions and to 
fix our notation. First of all we introduce the chirality, 
N5 = Nr — Nl, as the imbalance between right- and 
left-handed quarks. In QCD, change of chirality at zero 
as well as at finite temperature can be related directly 
to the topology of nonperturbative gluon configurations 
with a finite winding number, Qwi vi a the integrated 
Ward identity, 

N 5 = n 5 V = 2N f Q w , (1) 

where Nf is the number of flavors considered, and we 
have assumed that before the interaction with the gluon 
configuration, chirality was zero. In Equation ([1]) we have 
introduced the chiral density n§. The above equation is a 
consequence of the strong chiral anomaly; Qw is related 
to the topological charge of the given gluon configuration, 

At high temperature, a copious production of gluon 
configurations with nonvanishing winding number is ex- 
pected, see for example [4(| and references therein. These 
are the strong (that is, QCD) sphalerons. Therefore, be- 
cause of the existence of QCD sphalerons, chirality can 
be produced in the high temperature phase of QCD. 

The simplest way to treat quark matter with net chiral- 
ity in effective models, is to introduce an axial chemical 
potential, ^5, conjugated to chiral density [HI, E5I - I47I ]. 
At the lagrangian level, this amounts to add the chiral 
density operator, 

w^AV , (3) 

to the lagrangian density. This procedure is similar to 
the one usually adopted to study systems at finite baryon 
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density: in the latter case, we introduce a baryon chem- 
ical potential, /x, that induces a net baryon density, n. 
For example, in the very high temperatu re p hase of two 
flavor QCD, the following relation holds |4l| : 



n 5 



A 

3tt 2 



(4) 



which shows how a chiral chemical potential induces a 
net chirality in the system. In the general case, one needs 
to compute the relation between 725 and fi$ numerically, 
see [ill H(| and Section V of this article. 

Before going ahead, it is necessary to make few re- 
marks. Firstly, we are aware that the axial chemical 
potential cannot be considered as a true chemical po- 
tential. As a matter of fact, ^5 is conjugated to n$, or 
to the topological charge because of the Ward identity. 
The latter is a not conserved quantity in QCD because 
of the quantum anomaly. In the common picture of the 
QCD vacuum with ^5 = 0, the vanishing average value of 
the topological charge is understood as the formation of 
several local domains, each one characterized by a finite 
topological charge; the probability to create a domain 
with charge Qw is the same as the probability to cre- 
ate a domain with charge —Qw- As a consequence, in 
a volume much larger than the typical domain size, the 
average value of the topological charge sums up to zero. 
However, flucuations of the topological charge can change 
the value of the charge of any local domain. The formal 
role of the axial chemical potential is to reproduce the 
local domains of the QCD vacuum with a net topological 
charge as equilibrium states. Because of the fluctuations 
of the topological charge, this description is meaningful 
as long as the time scale is shorter than the inverse of the 
topological changing transition rate [4!| . This interesting 
theoretical question has been analyzed also in [52J , where 
it is claimed that a proper combination of the chiral den- 
sity and of a Chern-Simon term is a conserved quantity, 
and it is precisely this quantity that should be coupled 
to the axial chemical potential. We will devote a more 
detailed study of this problem to a next article. 

Another remark is that the theory with fi$ 7^ is a 
sign-free theory. It is well known that QCD with three 
colors suffers the sign problem: namely, the fermion de- 
terminant of QCD with three colors is complex at finite 
quark chemical potential, making the usual Montecarlo 
sampling of configurations in the Lattice simulations not 
possible when the quark chemical potential is larger than 
the temperature (see [|| for a review). On the other hand, 
the theory at finite [1$ does not suffer the sign problem. 
As a matter of fact, 75-0(1*5)75 = D'^^), where D cor- 
responds to the Dirac operator with ZX5 7^ 0. As a conse- 
quence, the fermion determinant is real and positive at 
fj>5 =/= 0, and grand canonical ensembles with finite ^ can 
be simulated on the Lattice [4l|. Indeed, some Lattice 
simulation at /Ltg 7^ has been already performed to study 
the chiral magnetic effect on the Lattice, see [H, H3|. In 
these references, a preliminary study of the phase struc- 
ture in the 1x5 — T plane is also addressed. 



III. THE MODEL 

In this Section, we describe the model that we use 
in calculations, namely the Nambu-Jona-Lasinio model 
improved with the Polyakov loop (PNJL in the follow- 
ing) [ll|. In the PNJL model, quark propagation in the 
medium is described by the following lagrangian density: 



C = q(i^D fl -m)q + G {qqf 



ini^rq)' 



(5) 



here q is the quark Dirac spinor in the fundamental repre- 
sentation of the flavor SU(2) and the color group; r cor- 
respond to the Pauli matrices in flavor space. A sum over 
color and flavor is understood. The covariant derivative 
embeds the QCD coupling with the background gluon 
field which is related to the Polyakov loop, see below. 

The PNJL model has access to the expectation value 
of the Polyakov loop, L, which is sensible to confinement 
or deconfinement properties of a given phase. In order 
to compute (L) we introduce a static, homogeneous and 
Euclidean background temporal gluon field, Aq = iA^ = 
i\ a A1, coupled minimally to the quarks via the QCD 
covariant derivative. Then W = Tr c exp (ifiA a A%) and 
L = W/3, where (3 = 1/T. In the Polyakov gauge, which 
is convenient for this study, Aq = iX^+iXg,^ 8 ; moreover, 
for simplicity we take L = from the beginning as 
in which implies A\ = 0. 

In our computation we follow the idea implemented 
in [2ll | . which brings to a Polyakov- loop-dependent cou- 
pling constant: 



G = g [l - aiLL 1 - a 2 {L 3 + (L f ) 3 )] 



(6) 



The ansatz in the above equation was inspired by [9( in 
which it was shown explicitly that the NJL vertex can be 
derived in the infrared limit of QCD, it has a non-local 
structure, and it acquires a non-trivial dependence on the 
phase of the Polyakov loop. This idea has been analyzed 
recently in [53|, where the effect of the confinement or- 
der parameter on the four-fermion interactions and their 
renormalization-group fixed-point structure has been in- 
vestigated. The numerical values of a\ and a 2 have been 
fixed in [2l[ by a best fit of the available Lattice data 
at zero and imaginary chemical potential of Ref . [5J, |55| . 
In particular, the fitted data are the critical temperature 
at zero chemical potential, and the dependence of the 
Roberge- Weiss endpoint on the bare quark mass. The 
best fit procedure leads to a.\ — a 2 = a = 0.2, which are 
the values we adopt in this article. 

In the one-loop approximation, the effective potential 
of this model is given by [47| 



d 3 p 



V=U{L,L\T) + a 2 G-N c N f Y J [P 



JV, fdhp_ { 

p hJ ( 27r ) 3 



(7) 
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where a = qq is a scalar collective field representing the 
mean field and the quantum fluctuations of the operator 
qq, and G is defined in Eq. ([6|); moreover, 



= \ K\p\s- 



(8) 



corresponds to the pole of the quark propagator, and 
m q = m — 2aG is the constituent quark mass; the index 
s denotes the helicity projection. In Eq. (|7|) we have 
introduced the functions 

T_ = 1 + 3Le~^ u '-^ + 3Z,t e -2«w.-/*) + e -W(u.-n) , 

(9) 

f+ = i + 3£t e -f + iLe -m^s+n) + e -m»s+») t 

(10) 

which are responsible for the statistical confinement 
property of the model at low temperature [ll[ . 

In right hand side of the first line of Equation ([7]) the 
momentum integral corresponds to the vacuum quark 
fluctuations contribution to the thermodynamic poten- 
tial. We treat the divergence in this term phenomeno- 
logically, introducing a momentum cutoff, A, in the vac- 
uum term; the numerical value of A will be then fixed 
by requiring that the quark condensate, the pion mass 
and the pion decay constant computed in the model are 
in agreement with the phenomenological values. Before 
going ahead, it is useful to remind that within the quark- 
meson model a renormalization procedure is feasible [l7| , 
the model itself being renormalizable. The renormaliz- 
ability of the model might be useful to remove the cutoff 
effects that, instead, appear in the NJL model (see the 
discussion in the next Sections). On the other hand, the 
renormalization program of the energy density at finite 
/i5 is not trivial, since a nonvanishing /15 induces further 
divergences in the theory, as it can be argued by an in- 
spection of the ^5— dependence of the vacuum energy, 
see for example Eq. (1141) . Hence, proper renormalization 
conditions should be adopted. 

The potential term U in Eq. (J7J is built by hand in 
order to reproduce the pure gluonic lattice data with 
N c = 3 [Hj]. We adopt the following logarithmic form, 



U[L,L,T] = T 4 <^ - 



a(T) 



LL 



+ b(T) ln[l - 6LL + 4(T 3 + L 3 ) - 3(LL) 2 ] j , 

(11) 

with three model parameters (one of four is constrained 
by the Stefan-Boltzmann limit), 

(12) 



'i(T) = n lt - ,i L ( +a 2 



i ( y 



(13) 



The standard choice of the parameters reads ao = 3.51, 
ai = —2.47, 0,2 = 15.2 and 63 = —1.75. The parameter 



To in Eq. (fTTj) sets the deconfinement scale in the pure 
gauge theory. In absence of dynamical fermions one has 
To = 270 MeV. However, dynamical fermions induce a 
dependence of this parameter on the number of active 
flavors [16}. For the case of two light flavors to which 
we are interested here, we take To = 190 MeV as in [2lj . 
Also for the remaining parameters we follow [2l[ and take 
A = 631.5 MeV, m = 5.5 MeV and g = 5.498 x 10~ 6 
MeV' 2 . 

We notice that the PNJL model considered here, which 
is dubbed Extended-PNJL in [2l|, has been tuned in or- 
der to reproduce quantitatively the Lattice QCD thermo- 
dynamics at zero and imaginary quark chemical poten- 
tial. Hence, it represents a faithful description of QCD, 
in terms of collective degrees of freedom related to chiral 
symmetry breaking and deconfinement. 



IV. PERTURB ATIVE ANALYSIS 

In the case of small (15, we can make some analyti- 
cal and semianalytical estimate of the effect of the ax- 
ial chemical potential on quark condensation, and on the 
critical line. In this Section we restrict our analysis to the 
pure NJL model case, corresponding to take L = = 1 
in Eq. ([7]), and ol\ — oti — in Eq. ([5]). This simpli- 
fies the numerical analysis. On the other hand, the re- 
sults obtained here will not be modified qualitatively by 
the Polyakov loop, since the effect of the latter is just 
a suppression of colored states below T c . Moreover, we 
work in the chiral limit; this simplification allows to de- 
fine rigorously the chiral phase transition, and compute 
unambiguously the critical temperature. 



A. Zero temperature: chiral condensate 

To begin with, we consider the zero temperature case, 
and we compute the shift of the chiral condensate in- 
duced by ^5, showing that the chiral chemical potential 
acts as a catalyzer of chiral symmetry breaking. In fact, 
the 115— dependent zero temperature correction to the ef- 
fective potential is given by 



_ N c Nf 2 2 p f m i\ 



(14) 



where we have assumed (i§ <C m q ; in Eq. (|14p we have 
introduced the function 



F(x) = loj 



l + VTT. 



1 



(15) 



It is easy to prove that F(x) is always positive, thus mak- 
ing V]° < 0. Hence, the energy density of a broken phase 
and /15 ^ is smaller than that of a phase with unbroken 
symmetry and the same value of the chiral condensate. 

We can use Eq. (jT4"f to analyze the perturbative solu- 
tion of the gap equation at zero temperature. As a matter 
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of fact, for small values of ^5 we can look for a solution 
of the gap equation in the form m q — fh 



5m q , where 



satisfies the gap equation at fi 5 = 0, and Sm q corre- 
sponds to the /is-dependent contribution. An elementary 
computation shows that 



Sm a 



J_ ( dV i\ 
Ml \dm q ) 



(16) 



where we have defined M„ 



(d 2 V /dml) m ^ and 

Vo = V(nz, = 0). By definition, M% > 0, since m q 
corresponds to a minimum of the effective potential at 
fj,5 = 0. Moreover, for m, C A we find (OV® / dm q ) oc 
m q log (m q /A) which is negative (we have verified that 
the sign does not change as long as m q < A). We con- 
clude that 5m q > 0, showing that /X5 favors the sponta- 
neous breaking of chiral symmetry at zero temperature. 
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FIG. 1. Function G of Eq. (|18|) versus temperature, for two 
values of the constituent quark mass, m q . Orange solid line 
corresponds to m q = 50 MeV; Indigo dashed line corresponds 



to the case m„ 



350 MeV. 



B. Finite temperature: chiral condensate 

We can extend the analysis of the previous subsection, 
to the case of finite temperature. For the perturbative 
solution of the gap equation, the derivative of the effec- 
tive potential with respect to the quark mass is needed, 
see Eq. (fT6|). At finite temperature, the /is-dependent 
contribution of the effective potential is Vi = V® + 
with Vf 1 given by Eq. (fLl}. Thus, Eq. ^ is replaced by 



Sm„ = — 



1 / dV? dV? 
M% \dm q dm q 



(17) 



The expression of is not informative, thus it is not 
necessary to report it here. Its derivative with respect to 
the quark mass is more interesting, 



dm„ 



N c N ft i 2 5 G(m q ,T) 



(18) 



In the above equation, G(m q ,T) corresponds to a func- 
tion defined in terms of a convergent numerical integral. 
In Fig. [TJ we plot G as a function of temperature, for two 
different values of the quark mass. We find that G > 0. 
As a consequence, because of Eq. (|T7l) , the /^-dependent 
contribution of the thermal fluctuations tends to reduce 
the value of the constituent quark mass (i.e., of the chiral 
condensate). 



Critical line 



We compute in this subsection the critical tempera- 
ture, T c , as a function of /Z5, and show at the same time 
that the chiral phase transition is of the second order. 
To this end we perform an expansion of the effective po- 
tential near T c . In fact, the order parameter around T c 
is small enough that a series expansion in powers of a/T 
is justified. Since we restrict ourselves to the case of 



the pure NJL model, there is no dependence from the 
Polyakov loop in the NJL coupling constant, G. 
One can write 



V = V 



Oti 4 



(19) 



where the dots denote higher order terms; Vq is the po- 
tential at a — 0: it is independent on the condensate, 
thus it is just a number which does not affect the physics 
of the problem. We notice that because of our defini- 
tion of the a field, which has the dimension of a cubic 
mass, we have extracted the proper powers of the NJL 
coupling constant from the definitions of the coefficients, 
in order to give the common mass dimension to the coef- 
ficients themselves. At the second order transition point 
T = T c one has 012 = and 04 > 0. Therefore, to de- 
termine the critical temperature as a function of /i 5 it is 
enough to determine the zeros of 0,2 m the ^5 — T plane. 
Next, we compute q?4 to check that the transition points 
correspond to a second order phase transition. 

Computing the second derivative of the effective po- 
tential in Equation ([7]) we find, after some algebra, 



OL2 = «2,0 + "2,2^5 



where 



«2,0 



«2,2 



2 
G 

4N r N 



2N C N 



log 



A 2 

T 
A 



2Ayv/ r 2 



(20) 

(21) 
(22) 



and the numerical factor c = 1.36 arises from a conver- 
gent numerical integral. 

At small /is the equation ct2 — 0, which determines the 
critical temperature, can be solved perturbatively: from 
the condition 0:2.0 = we find 



rpO _ 



3A 2 



N r N f G 



1/2 



(23) 
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using the parameters of the model we find T® = 174 MeV. 
Then, writing T c = T® + ST and solving 0:2.0 + a 2,2^g = 
for ST we find 



T B = T? ( 1 



3/*§ 



4N c N f (T°) 



r«2,2 



(T c °) 



(24) 



Using the parameters given above we find 0:2,2(^0') = 
0.17, which implies that the critical temperature is a de- 
creasing function of Hence, a finite /is favors the 
disordered phase at finite temperature. 

The last step is to check the sign of 04, to be confident 
that the transition is of the second order. Since we make 
a perturbative analysis at /15 -C T, we may assume /ig 
infinitesimal of the same order of a, hence terms of the 
order of C(^|cr 4 ) can be regarded as 0(a 6 ) and can be 
neglected. At the leading order, it is therefore enough to 
compute the coefficient at /ig = 0, 



A8N c N f / A 

«4 = ^^ bg- 



(25) 



where b — 0.90 arises from a convergent integral com- 
puted numerically. Using the values of A and Tj? we find 
that a,\ — 11.35 at T — T°; being it positive, the phase 
transition is of the second order. 

Before going ahead, it is useful to summarize the re- 
sults of this Section. Eq. (|17p represents the correction to 
the solution of the gap equation due to /ig 7^ 0. The first 
and second addenda in the right hand side of the equation 
correspond to the zero temperature and finite tempera- 
ture contributions, respectively. At zero temperature, 
only the first contribution survives, and it is negative, 
leading to a positive shift of the constituent quark mass. 
On the other hand, the contribution of thermal fluctua- 
tions is positive. Thus at a given temperature, there is a 
competition between the effect of /ig on the vacumm and 
thermal contributions to the gap equation. In passing, 
we notice that if vacuum fluctuations were neglected in 
the thermodynamic potential, then the thermal contribu- 
tion in Eq. (1171) would lower the constituent quark mass 
at finite temperature as /is is switched on. Finally, the 
results in Equations (|2"21 , (j2~4")l and (|2"5j) show that at zero 
and small /ig, the chiral phase transition is of the second 
order, and the critical temperature is a decreasing func- 
tion of the chiral chemical potential. This scenario for 
the critical line is in agreement with the results of (45l — 
H3|, where a numerical procedure of minimization of the 
one-loop potential has been adopted. 



V. GRANDCANONICAL AND CANONICAL 
PHASE DIAGRAMS 

In this Section, we describe the expected phase struc- 
ture of QCD at finite temperature and axial chemical 
potential, as predicted by the PNJL model. The grand- 
canonical phase diagram, already discussed in [47] . cor- 
responds to the map of the different phases in the /ig — T 



plane. On the other hand, for the case of the canoni- 
cal ensemble phase diagram, /ig is replaced by the chiral 
density n 5 , the latter being defined as 



n 5 



(26) 



where the derivative has to be computed at the global 
minimum of the grand potential. 

The computation of 715 in the model needs some care. 
In particular, at m q 7^ the divergence of the vacuum 
energy at finite /is is transmitted to rig. As a matter of 
fact, 715 can be formally split into a vacuum and a thermal 
parts. The thermal part is convergent, and its derivative 
is finite as well, therefore it does not lead to any difficulty. 
On the other hand, the T = contribution is divergent, 
and its derivative with respect to //g is divergent as well 
in the case m q 7^ 0. This can be realized easily if we take 
for a moment the limit A»/is, m qi in which we find 



N r N 



n 5 



2vr 2 



2 3 



, 2A 
log 1 

m„ 



(27) 



In the case m q = we obtain the result of [4l|, which is 
not affected by the ultarviolet divergence. Because of the 
latter, we cannot take the limit A — ¥ 00 in the calculation 
of the phase space integrals which are involved in the ex- 
pressions of 7ig. Therefore in the numerical computation 
of ri5 , we cutoff the vacuum contribution at p = A unless 
otherwise stated, for internal consistency. 

Before the discussion of the results, it is useful to com- 
ment briefly on the expected effect of the finite cutoff 
on ng. If /is 3> mo, which is appropriate in the quark- 
gluon plasma phase, we can neglect the mass term in the 
vacuum contribution of the fermion determinant to rig, 
which we call Co- 



Co 



N c N f 
3tt 2 



[/^(A-^+A^g-A)] 



(28) 



We notice that the term proportional to A 3 in the above 
equation does not appear in Eq. (|27p since in the latter 
we have taken the limit A > /J5. As soon as /ig > A, 
the vacuum contribution to the chiral density saturates 
because of the existence of a cutoff in the theory, as shown 
by Eq. (|28[). In the numerical calculations, we expect 
that the saturation effect will be softened by the thermal 
fluctuation contributions. As a consequence we expect 
to measure some mild saturation effect as /ig « A. This 
saturation is observed on the Lattice as well, see [5(| ■ 
In last analysis, this is not a serious trouble, since we 
expect the model to be consistent only when masses are 
smaller than the cutoff, in this case for regime /ig < A. 

In Fig. [5] we plot the normalized chiral density, pw = 
n$/2Nf, as a function of /is for several values of the tem- 
perature. In the Figure, T c — 173.9 MeV corresponds 
to the critical temperature for deconfinement and chiral 
symmetry restoration at /Z5 = 0. At large values of /is we 
measure the expected saturation of pw ■ This saturation 
is not physical but a mere artifact of our regularization 
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100 200 300 400 500 600 700 

fi 5 (MeV) 

FIG. 2. Normalized chiral density, pw = n$/2Nf, as a func- 
tion of the chiral chemical potential, at several values of the 
temperature. From right to left, temperature is equal to 
0.4T C , 0.6T C , 0.8T C , 0.9T C , T c5 , 1.1T C respectively. 



scheme. Hence we will consider trustable only the results 
in the range p,^ < A; when the results outside the above 
specified range are shown, it is done only for complete- 
ness. 

Next we turn to the phase structure. The grandcanon- 
ical phase diagram has been computed in [43]. It is use- 
ful, however, to briefly summarize here the results of [47| . 
in order to facilitate the comparison with the canonical 
phase diagram. In the upper panel of Fig [3] we plot the 
phase diagram in the p.§—T plane. The indigo dashed line 
corresponds to the chiral and deconfinement crossover; 
the solid line denotes the first order phase transition. The 
indigo dot is the critical endpoint. The scale T c = 173.9 
MeV corresponds to the critical temperature at [i$ = 0. 

The critical line in the phase diagram in Fig. [3] is iden- 
tified with the peak of the derivative dL/dT. Within few 
MeV, we have found that the latter coincides with the 
location of the peak of \da/dT\, for the values of pi$ an- 
alyzed in this work. Thus, within the model ad hand, 
the deconfinement and the chiral symmetry restoration 
take place simultaneously. In last analysis, this is related 
to the explicit dependence of the NJL coupling on the 
Polyakov loop, and to the mechanism that leads to de- 
confinement in the model. As a matter of fact, at finite 
temperature the one- and two-quark contributions to the 
grand potential act as efficient Z3— breaking terms, which 
in turn favor the L ^ state. As a consequence, the ef- 
fective value of the NJL coupling which is responsible 
for the spontaneous chiral symmetry breaking is reduced 
with respect to the zero temperature case, as soon as 
L ^ 0, see Eq. ([6]). Thus, it is expected that deconfine- 
ment in the extended PNJL model implies chiral sym- 
metry restoration, as it is confirmed by the numerical 
results. The inverse assertion, namely that chiral sym- 
metry restoration implies deconfinement, is not true and 
indeed it is not necessarily realized in the model, unless 
an explicit backreaction on the Polyakov l oop effective 
potential is introduced, see for example [l4], [HI, HH ■ 
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FIG. 3. (color online). Upper panel: phase diagram in the 
grandcanonical ensemble. The indigo dashed line corresponds 
to the chiral and deconfinement crossover; the solid line de- 
notes the first order phase transition. The indigo dot is the 
critical endpoint. The scale T c = 173.9 MeV corresponds to 
the critical temperature at p$ = 0. Lower panel: phase di- 
agram in the temperature-chiral density plane. Solid lines 
correspond to the boundaries of the mixed phase. On the left 
line the volume fraction, c, of hadron phase is c = 1; on the 
right line we find c = 0. Dashed line denotes the chiral and 
deconfinement crossovers. Dot-dashed line corresponds to the 
values of T(pw) at which c = 1/2. In the portion of the phase 
diagram below the dotted line we find /15 > A. Finally, the 
indigo dot denotes the critical endpoint. 



The grandcanonical phase diagram can be translated 
to a canonical one, by replacing /15 with pw- This pro- 
gram is easily accomplished once we use Eq. (1261) to com- 
pute pw, once the physical values of (qq) and L are 
known. As already explained in the Introduction, the 
use of the term canonical has to be taken as an abuse 
of nomenclature. Indeed n§ is not a conserved quantity, 
because of the fluctuations of the topological charge in 
QCD and the Ward identity. Thus, discussing about 715 
is meaningful only if the time of observation of the sys- 
tem is smaller than the typical time needed to observe 
a fluctuation of the topological charge, the latter being 
related to the inverse of the sphaleron transition rate. 

In the lower panel of Fig. [3J we plot the phase diagram 
in the pw — T plane. In the figure, the indigo dashed 
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line denotes the chiral and deconfinement crossovers (the 
two coincide within numerical uncertainties within this 
model calculation). The indigo dot denotes the critical 
endpoint. Its coordinates in the phase diagram are 

pi = 0.518 tuT 3 , T 5 C = 167 MeV . (29) 

The orange solid lines correspond to the boundaries of 
the mixed phase, which develops in correspondence of 
the first order phase transition line of the grancanoni- 
cal phase diagram. The composition of the mixed phase 
at the temperature T can be evaluated easily as follows. 
Given the total charge density pw, then the volume frac- 
tions of the ordered, c, and disordered, d, phases have to 
satisfy cp\ + dp-2 = pw with c + d = 1; here p\^ denote 
the values of the topological charge density on the two 
boundaries of the mixed phase region. It follows then 



P2 - Pi 

on the left line the volume fraction c = 1, corresponding 
to a homogeneous pure ordered phase; on the right line 
c = 0, corresponding to a pure disordered phase. For 
reference, in Fig. [3] we plot a pink dot-dashed line which 
corresponds to the values of T(pw) at which c = 1/2. 
Finally, in the portion of the phase diagram below the 
green dotted line we find ^5 > A, in which cutoff artifacts 
are relevant. At small temperatures the green line is 
almost vertical, because the topological charge density 
turns out to be very insensitive of temperature, see Fig. [2] 

VI. INFLUENCE OF THE QUARK MASS AND 
OF THE VECTOR INTERACTION 

In this Section, we wish to investigate the quantitative 
effect of a nonphysical bare quark mass, as well as of the 
presence of the vector interaction, on the critical line in 
the p^ — T plane. In particular, we wish to compute the 
evolution of the critical endpoint coordinates as a func- 
tion of the bare quark mass, and of the coupling strength 
in the vector channel. 



A. The effect of the bare quark mass 
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FIG. 4. (Color online). Evolution of the critical endpoint 
coordinates as a function of the bare quark mass. Solid or- 
ange line corresponds to the critical value of the axial chem- 
ical potential; dashed indigo line corresponds to the critical 
endpoint temperature. Axial chemical potential and temper- 
ature are measured in units of the endpoint coordinates at 
m = m ,ph ys - 



orange line corresponds to the critical value of the axial 
chemical potential; the dashed indigo line corresponds to 
the critical endpoint temperature. Axial chemical poten- 
tial and temperature are measured in units of the end- 
point coordinates at mo = mo.phys with mo,phys = 5.5 
MeV, which corresponds to ps 139 MeV. For compar- 
ison, at the largest value of quark mass considered here 
we find run m 400 MeV. 

Our results show that the critical endpoint tempera- 
ture, Tg , is not so much affected by the quark mass. As 
we will show in the next Section, this inertia of T5 is also 
present when we switch on a vector interaction. On the 
other hand, the critical endpoint axial chemical potential, 
/ig, is strongly affected by the quark mass. In more detail, 
the larger value of mo, the larger p,\. This is quite easy 
to understand naively, since the bare quark mass turns 
the chiral transition to a crossover. At p 5 = 0, the chiral 
crossover is smoothed as mo is increased. Therefore, it is 
natural to expect that a larger value of ^5 is necessary to 
experience a first order phase transition, if m > m .phy S - 
This might partly explain why the critical endpoint CP5 
is not yet detected in Lattice simulations [49|, [5(| • 



In this Section, we study the effects of the bare quark 
mass on the location of the critical endpoint CP5. Esti- 
mation of the effect of the bare quark mass is very im- 
portant, because Lattice simulations at /its 7^ nowadays 
are performed with a numerical value of the quark mass 
which is larger than the physical value. In Lattice simula- 
tions, the non physical value of the quark mass is reflected 
into a non physical value of the pion mass. In [4!| [o0| . 
the value of the pion mass in the vacuum is m w ss 400 
MeV. 

In Fig.@]we plot the critical endpoint coordinates as a 
function of the bare quark mass. In the figure, the solid 



B. The role of a vector interaction 

In this Section we briefly comment on the role of a 
vector interaction on the phase structure of the model. 
To this end, we add to the lagrangian density the term 



(^V) 2 + (^75^) 2 



(31) 



we do not include the interaction in the triplet channel, 
since at the one-loop level it gives rise to terms which 
couple the isospin density, 5n = rid — n u , to the isospin 
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density operator; these are not relevant as long as we do 
not introduce the isospin chemical potential. 
At the mean field level one has 



C v = Gyn 2 - 2GVm/W 

+G v nl - 2G v n 5 j} 1 ° 1 5 iP ; 



(32) 



as usual, we have defined n = (ifi'ip) and = (V>'7 5 'i/ ; )- 
Thus at the mean field level, the chemical potentials for 
quark number density and for chiral density are shifted 
because of the vector interaction, 



jit — > /j, — 2Gyn , 
/i 5 fi 5 - 2G v n 5 . 



(33) 
(34) 



The previous equations are useful to grasp the effect 
of Gv 7^ on the phase structure. As a matter of fact, 
if Gy > then quarks propagate in a thermal bath at 
an affective chemical potential which is smaller than the 
true chemical potential, because of Eq. (|33|) . As a con- 
sequence, the transition at finite \x is smoothed, and the 
critical endpoint moves to higher values of the chemi- 
cal potential. On the same footing, because of Eq. AMI) , 
the effective chiral chemical potential is smaller than n§ . 
Therefore we expect that the critical endpoint coordi- 
nate moves to higher values of ^5 compared to the case 
G v = 0. 

This reasonings are confirmed by our concrete numer- 
ical computations. In the case of Gv 7^ 0, the value 
of n 5 has to be computed self-consistently by means of 
the number equation (|26[) , which has to be solved, at any 
given value of /Z5 and T, together with the gap equations, 



on _ dn 



(35) 



In Fig.[5]we plot the expectation value of the Polyakov 
loop in the pseudocritical range as a function of temper- 
ature, for several values of the strength of the coupling 
in the vector channel, at /15 = 300 McV. At Gy = 0, 
this value of is slightly above In the figure, solid 
line corresponds to Gy = 0; dashed line corresponds to 
Gv = 0.25G; dot-dashed line corresponds to Gv = 0.5G; 
dotted line corresponds to Gy — 0.75G; finally, short 
dashed line corresponds to Gy — G. As expected, the 
role of the vector interaction is to smooth the phase tran- 
sition in comparison with the latter at Gy = 0. In- 
deed, the first order phase transition at /X5 = 300 MeV 
at Gy — is turned by the vector interaction into a 
crossover at the same value of /X5. 

In Fig. [5] we plot the critical endpoint coordinates, 
and T c , as a function of the ratio Gy/G. Blue dashed 
line corresponds to the critical temperature; red solid 
line corresponds to the critical chiral chemical potential. 
We find that the temperature of the critical endpoint 
is poorly affected by the value of Gy in the range that 
we have examined; on the other hand, the critical value 
of the chiral chemical potential is quite sensitive to the 
vector coupling. 




162.5 165 167.5 170 172.5 175 177.5 180 
T (MeV) 

FIG. 5. {Color online). Expectation value of the Polyakov 
loop in the pseudocritical range as a function of temperature, 
for several values of the strength of the coupling in the vector 
channel, at /15 = 300 MeV. Black solid line corresponds to 
Gv = 0. Green dashed line corresponds to Gv = 0.25G. Blue 
dot-dashed line corresponds to Gv = 0.5G. Orange dotted 
line corresponds to Gv = 0.75G. Finally, indigo short dashed 
line corresponds to Gv = G. 
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FIG. 6. [Color online). Critical endpoint coordinates as a 
function of the ratio Gv / G. Indigo dashed line corresponds 
to the critical temperature; orange solid line corresponds to 
the critical chiral chemical potential. Axial chemical poten- 
tial and temperature are measured in units of the endpoint 
coordinates at mo = mo.Phys- 



It is useful to comment about the similarities between 
the phase structure at finite /x 5 , which turns out from our 
computation, and the one well established (within effec- 
tive models) at finite baryon chemical potential, /i. In 
this article, we have considered only one type of conden- 
sate, namely the chiral condensate, which characterizes 
the symmetry breaking pattern at finite T and /X5 in the 
fermion sector of our model. Besides a smaller curva- 
ture of the critical line in the case of the chiral chemical 
potential, compared with that at finite fi, we do not find 
qualitative difference in the phase structure. Also, the ef- 
fect of the vector interaction, as well as of a finite current 
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quark mass, is very similar, qualitatively, in the cases of 
finite \x and finite /15. Some difference between the two 
phase diagrams might arise from different types of con- 
densates. For example, an important point that we have 
not considered in our article, for the sake of simplicity, 
is the introduction of a diquark condensate, which might 
appear at finite mu5. For the case of QCD at finite /z, 
such a condensate is expected to be developed at very 
large /x [62l - [64j ]. However, this possibility deserves further 
study and needs to be checked numerically by dynamical 
computation of the condensates. We leave this point to 
a future project. 
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VII. TOPOLOGICAL SUSCEPTIBILITY 

In this Section we compute the topological suscepti- 
bility at finite temperature, in presence of a background 
chiral density. Topological susceptibility in QCD is de- 
fined as the correlator of the topological charge at zero 
momentum; it can be computed from the QCD partition 
function via the relation 



X 



d 2 n 



(36) 



1=0 



where the 0-angle is introduced by adding the following 
term to the lagrangian density: 



) 9 riivpo pa pa 

64tt 2 pa ' 



(37) 



It is well known that one can get rid of the term in 
Eq. (|3"T|) in the QCD action by means of a chiral ro- 
tation. However, after the chiral transformation, the 9 
dependence of the QCD partition function appears ex- 
plicitely in the quark part of the QCD action. A model 
study having some overlap with our study can be found 
in [ijj, where the full topological susceptibility in pres- 
ence of a back grou nd of topological charge is computed. 
In the case of [19| . the topological charge is introduced 
by adding the 9— term to the lagrangian density, which in 
turns acts as a source for the topological charge. The rel- 
evance of a finite 9 in [l|| lies on the possibility that for 
temperatures higher than Aqcd, sphaleron transitions 
induce a spacetime dependent 0-angle [5lj . 

A low energy relation connects the vacuum chiral con- 
densate to the topological susceptibilty; in the theory of 
Nf light flavors the relation reads [HD, [57} 



O 



A 



QCD 



(38) 



where (qq) is the common value of the quark condensates 
for the light quarks. Equation (|38[) shows that the topo- 
logical susceptibility is proportional to the product of the 
light flavor masses; thus, in the equation it is manifest 
the fact that in a theory in which at least one massless 



FIG. 7. Topological susceptibility as a function of temper- 
ature, for several values of /is. Black data correspond to 
fis = 0; orange data correspond to lis = 300 MeV. The green 
dotted lines correspond to the right hand side of Eq. 



flavor exists, one has \ = 0- This was proved also by 
Ward identities in (58j . 

As stated before, a chiral rotation transmits the 9- 
dependence of the QCD action to the quark sector. After 
the rotation, the quarks acquire a complex mass term. 
Besides, the condensates that are not invariant under 
the axial rotation are mixed among themselves. This 
scenario can be implemented within the effective mod- 
els to compute 9— dependent quantities, see for exam- 
ple [Hll . Beside this, the 0-dependent action can be used 
to compute \ by virtue of Eq. (|3"o) . A detailed discussion 
about the implementation of the 9— dependent lagrangian 
within the chiral model with the Polyakov loop can be 
found in [l9|, [20|. In this article it is enough to men- 
tion that, in order to introduce the 9 angle in the model 
and compute the topological susceptibility, it is enough 
to change M — tuq — 2aG into the grand potential with 



in,, s i's [ - 1 — 2<jG 



+ to 2 . sin 2 I - 



(39) 



It is worth to mention here that generally speaking, 
at finite 9 other condensates might develop (and indeed 
they do develop, see [III H(|). The presence of these 
condensates makes the substitution (|3"9"|) not sufficient for 
the complete treatment of the 9 angle within the chiral 
models. However, such new condensates vanish in the 
theory at 9 = 0, to which we are interested; neglecting 
them from the beginning does not change the derivative 
of f2 at 9 = 0. Therefore, the replacement in Eq. (f39|) 
is sufficient for the purpose of computing the topological 
susceptibility within this model. 

In Fig. [7] we plot the topological susceptibility as a 
function of temperature, for several values of /X5 . At zero 
temperature and chiral chemical potential we find 



X = (79.97 MeV) 4 



(40) 
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which is in agreement with the large N c prediction in 
Eq. ([55]). which gives x = (80.3 MeV) 4 with two light 
flavors at T = /is = 0. 

The relation (|38p is actually satisfied within the chiral 
model also at T — and /is ^ 0. Indeed, a straightfor- 
ward computation shows that 



where 



1=G 



m N c Nf 
' 2 2tt 2 



p 2 dp 



al 



(41) 



V(p + /^) 2 + (™o-2ctG) 



A*5 



(42) 

moreover, from the conditions dfl/da — and a = 2(uu) 
we have 



(««) 



^^ jX I Q ( m ° 
2tt 2 \ m„ 



(43) 



Here m q corresponds to the constituent quark mass. A 
comparison between Eqs. (|41 [) and ([43"[) leads to 



(44) 



m °/- \ 
X = — 2~{uu) 



. m„ 



which is in agreement with Eq. (|38l) when the latter is 
computed for two degenerate flavors. On the other hand, 
at finite temperature we measure some discrepancy be- 
tween the theoretical prediction given by Eq. ([38]) . cor- 
responding to the green dotted line in Fig. [7j and the 
numerical results obtained within the PNJL model (solid 
lines in the same figure). 



Secondly, we have discussed the phase diagram in the 
canonical ensemble formulation, in which /is is replaced 
by the chiral density, n^. This might be interesting be- 
cause ns in QCD is connected locally to the topological 
charge density, because of the integrated Ward identity. 
As a consequence, it might be of interest to predict the 
numerical value of the topological charge density along 
the critical line, as well as at the critical endpoint. 

Thirdly, we have computed the effect of the vector in- 
teraction, as well as of the finite current quark mass, on 
the location of the critical endpoint. Such aspects should 
be taken into account if a comparison with the Lattice 
data is desirable, and extend the study started in [47 1. 

As a final investigation, we have computed the full 
topological susceptibility (i.e., which takes into account 
both the pure gauge and the dynamical fermion contri- 
bution) at zero as well as at finite /is, as a function of 
temperature. We find that the Di Vecchia-Leutwyler- 
Smilga-Veneziano (DLSV) relation [56|, [57| is satisfied at 
finite /is in the confinement phase of the model. At large 
temperature, above the critical temperature, we measure 
a deviation from the DLSV relation, both at zero and 
at finite /is; this can be understood within the model, 
since terms of the order of mo / m q which are negligible in 
the confinement phase, become important in the quark- 
gluon-plasma phase. 

It is interesting to ask wether the work presented here 
can be improved. In our opinion, several directions are 
possible for future research. As a first step, it would be 
important to investigate from a theoretical point of view, 
how to couple correctly /is to a conserved quantity in 
QCD. Along this line, the work in Ref. [52j seems quite 
enlightening. Secondly, it is of interest to include the 
possibility of other condensates, among them a diquark 
condensate, in analogy to the situation of QCD at very 
large density. We plan to report on these, as well as 
related, topics in the next future. 



VIII. CONCLUSIONS 
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